Tuning laser-induced bandgaps in graphene 
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Could a laser field lead to the much sought-after tunable bandgaps in graphene? By using Floquet 
theory combined with Green's functions techniques, we predict that a laser field in the mid-infrared 
range can produce observable bandgaps in the electronic structure of graphene. Furthermore, we 
show how they can be tuned by using the laser polarization. Our results could serve as a guidance 
to design opto-electronic nano-devices. 

PACS numbers: 73.23.-b, 72.10.-d, 73.63.-b 
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More than a century ago, the use of alternating cur- 
rents (ac) sparked a revolution that changed our mod- 
ern world. Today, the use of ac fields has reached the 
nanoscale [H, 0. Here, the interplay between the quan- 
tum coherence of the electrons, inelastic effects and dy- 
namical symmetry breaking offers fascinating opportu- 
nities for both basic research and applications. The co- 
herent destruction of tunneling Q and quantum charge 
pumping 04^5 the generation of a dc current at 
zero bias voltage due to quantum interference Q, are 
prominent examples of the wealth of phenomena driv- 
ing this rapidly advancing area of research. By means 
of time-dependent gatings 0, Q , surface acoustic waves 
flfjj or illumination with a laser field [Til ], experiments 
have probed different facets of ac transport and light- 
matter interaction. On the other hand, theoretical in- 
sights keep opening excitin g di rections, from interaction- 
induced quantum pumping jl2| to everlasting oscillations 

Among many other benefits, the advent of graphene 
[T3] and carbon nanotubes 15 1 provide an outstanding 
arena for the study of ac transport in highly coherent, 
low dimensional systems. Adiabatic quantum pumps 
[Ttj \Tf\ . ac controlled Fabry-Perot resonators [TH] and 
photodetectors are among the new breed of carbon- 
based devices. Furthermore, a few remarkable studies on 
Dirac fermions interacting with linearly [2(3, H3 and cir- 
cularly [22|, [23[ polarized monochromatic light pointed 
out striking non-perturbative non-adiabatic effects: a 
laser field could lead to the opening of dynamical gaps 
pfjl ] in graphene. For circularly polarized light, a further 
gap was shown to develop at the Dirac point |22j . Many 
open questions remain: which experimental setup would 
unveil these phenomena and how could we tune them? 

In this Letter, we aim at ellucidating these questions 
by analyzing the interaction between electrons and a 
monochromatic laser field of arbitrary polarization in 
graphene. Using a Floquet approach one finds that a 
laser of frequency Q induces the lifting of degeneracies be- 
tween (Floquet) states of the combined electron-photon 



system. Here, a careful analysis allows for the tuning of 
a feasible parameter range (laser frequency, power and 
polarization). Our predictions show that these effects 
are within the reach of mid-infrared laser technology in 
a transport setup, thereby opening promising prospects 
for graphene-based opto-electronic devices. 

In graphene, the low energy electronic states contribut- 
ing to transport are close to the Dirac points K and 
K'. Since we consider a clean sample and given that 
the ac field does not introduce any inter-valley coupling, 
we can describe both points separately. In the K-valley, 
those states can be accurately described by the k • p ap- 
proximation through the envelope wave-function = 
(^A, ^b) T , where the two components refer to the inter- 
penetrating sublattices A and B [24[. The time periodic 
electromagnetic field, with period T = 2n/il, is included 
as a monochromatic plane wave traveling along the z- 
axis, perpendicular to the plane defined by the graphene 
sheet. The magnetic vector potential is thus written as 
A(t) = Re{A e~ in *}, where A = A (l,e iv ) refers to 
the intensity and polarization ip of the field. For this 
choice of parameters, ip — yields a linearly polarized 
field A(t) = Aq cos fit (x + y) while <p = ir/2 results in a 
circularly polarized field A(i) = Ao(cosf2t x + sin ilt y). 
Therefore, the graphene electronic states in the presence 
of the ac field are encoded in the Hamiltonian: 



H(t) = v F a ■ [p - eA(t)} . 



(1) 



Here vp ~ 10 6 m/s denotes the Fermi velocity and 
o = {?x 7 &y ) the Pauli matrices describing the pseu- 
dospin degree of freedom. The suitable wavelength of 
the laser field lies within the mid-infrared region, i.e. 
A ~ 9 fim, HQ ~ 140 meV and we choose two intensity 
values I ~ 32 mW/ /im 2 and I ~ 130 mW/ /1m 2 such that 
eAoVF/tiQ — 0.17 and 0.33 respectively. 

As will be made clear later on, a correct description 
of our problem crucially requires a solution valid beyond 
the adiabatic approximation. Hence, the Floquet the- 
ory Q represents an appropriate approach to describe 



such electron-photon scattering processes. The resulting 
picture is that of an effective time-independent Hamil- 
tonian in a higher-dimensional space, the so-called Flo- 
quet space, defined as the direct product between the 
usual Hilbert space and the space of 27r/f2-time peri- 
odic functions. This space is spanned by the states 
{|k, n)_j_}, where k is the electronic wave vector, n is the 
Fourier index, and the subindex refers to the alignment 
of the pseudo-spin with respect to the momentum. In 
this basis, the problem to solve is identical to the time- 
independent Schrodinger equation with Floquet Hamilto- 
nian Hp = H—ihdt . The power of recursive Green's func- 
tion techniques [25| can be exploited to obtain both the 
dc component of the conductance and density of states 
from the so-called Floquet Green's functions [261 ]. 

Previous studies [2(| [HJ have considered lasers either 
in the far infra-red (M7 = 29 meV) [2(| or in the vis- 
ible range [22^ . In the first case, the predicted gaps 
where of about 6 meV and the authors considered the 
photocurrents generated in a p-n junction. On the other 
hand, for a laser in the visible range, the photon en- 
ergy of about 2eV is much larger than the typical op- 
tical phonon energies (170meV), severe corrections to 
the transport properties due to dissipation of the excess 
energy via clcctron-(optical)phonon interactions are ex- 
pected. Besides, appreciable effects in this last case re- 
quired a power above 1 W//im 2 , which could compromise 
the material stability. To overcome both limitations we 
quantitatively explore the interaction with a laser in the 
mid-infrared range (A = 5 — 10/im), where photon ener- 
gies can be made smaller than the typical optical phonon 
energy while keeping a much lower laser power. Fur- 
thermore, we show that the polarization, whose role was 
overlooked, can be used as a control variable. 

Figure [1] shows how the dc density of states (DOS) 
evolves as the polarization changes. Already for linear 
polarization (upper panel) one can find some surprises 
not reported before: i) close to the Dirac point, the DOS 
is enhanced as compared to the case without laser (red 
dotted line). The DOS increases linearly but with a dif- 
ferent slope, ii) Close to ±M2/2, stronger effects lead to 
a depletion of the states but without reaching a full gap. 
The counterpart of this depiction is given by the peaks 
in the DOS on each side of the depletion area. More 
accurate calculations based on a tight-binding model in 
a /im sized sample interacting with a linelarly polarized 
field confirm these findings (top panel). 

The successive panels in Fig. [T]show how these features 
change with the laser polarization: i) close to the Dirac 
point a gap opens and the structure of the side peaks is 
severely modified, ii) The depletion areas around ±ft£l/2 
transform into gaps whose widths are maximum for cir- 
cular polarization. A closer look at these figures show 
that the mentioned gaps are areas with a reduced but 
non-vanishing DOS. 

To rationalize the behavior observed in Fig. [T] we re- 
sort to the Floquet picture explained above. A scheme 
with the relevant states close to the Dirac point is shown 
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FIG. 1: (Color online) Scheme of the considered setup where 
a laser field with Kl = 140 meV is applied perpendicular to 
a graphene mono-layer. Density of states for (a) linear, (b) 
tp = 0.125-7T, (c) Lp = 0.3757T, and (d) circular polarizations. 
The black solid line is for I = 32mW/pm 2 while the red 
dashed line corresponds to I = 130mW/^m 2 . For linear po- 
larization these results are compared with those of a tight- 
binding calculation for a system with 5 x 10 4 channels (blue 
and purple lines). For reference the zero- field DOS is shown 
in solid grey. 



in Fig. [2j On the left panel, one can see the dispersion 
relation for the states {|k, n) ± } for n = (black solid 
line) and n = ±1 (red line), the scheme on the right 
panel shows the states represented by circles and their 
corresponding interactions with full lines. The effects of 
the ac field are expected to be stronger at the crossing 
points between these lines (whenever they have a non 
vanishing Hamiltonian matrix element), leading to the 
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FIG. 2: (Color online) (a) Scheme of the quasienergies as 
a function of k. The relevant crossing regions, marked by 
blue and red circles, yield the dynamical and central gap 
respectively, (b) Representation of the Floquet Hamilto- 
nian for the k • p approach. Here, the direct hopping vi = 
eAqv f ^ cqs a _|_ g-'VgjjjQ,) se ^ s the transition amplitude be- 
tween Floquet states conserving pseudospin. In contrast, the 
off-diagonal term vi = i eA <^" F (e~ lv cos a — sin a) introduces 
an inelastic back-scattering process that simultaneously en- 
ables both m and pseudospin transitions. 

opening of energy gaps at those points. From a geomet- 
rical argument one can see that the crossing of the states 
differing in one photon lies exactly at ±M2/2. These de- 
generacies are lifted by the ac field leading to the gaps 
observed in Fig. [T] This mechanism is reminiscent of an 
inelastic Bragg reflection. An estimation of the resulting 
gap gives: 

&k=n/2v F - eA v F y/l -cos(yj) sin(2a), (2) 

where a = tan -1 (k y /k x ). Equation [2] shows that the 
effect is linear in the field strength. Interestingly, after 
summing up over all the directions in the 2d k-space, we 
see that no net gap opens in the linearly polarized case 
since in the orientation a = ir/A both states become de- 
generate. However, as can be seen in Fig. [TJ there is a 
strong modification in the DOS that resembles the usual 
Dirac point for a Fermi energy around Ml/2. Chang- 
ing the polarization away from the linear case, one sees 
that a dynamical gap [20| opens and reaches its max- 
imum A max ~ 23meV (46meV) for I = 32mW/^m 2 
(130mW/rmi 2 ) in the circularly polarized case. 

Now, let us account for the phe nomenon occuring at 
the Dirac point. Oka and Aoki [22| predicted that circu- 
larly polarized light would induce a further gap around 
the Dirac point, in agreement with our Fig. QJd. Here we 
see that these strong modifications extend all the way up 
to linear polarization. A careful analysis shows that it 
is a higher-order effect. The leading contribution comes 
from the states connecting the degenerate ones at the 
Dirac point (there are four of these paths as can be seen 
by analyzing Fig. [2]-b). The calculated gap reads: 

A fc=0 ~ — Re{i;i?; 2 } = 2 — sirup, (3) 
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FIG. 3: (Color online) dc conductance through a graphene 
stripe of \\im times lfim in the presence of a linearly polarized 
laser as a function of the Fermi energy. The black solid line 
is for / = 32 mW / /im 2 while the red dashed line corresponds 
to 7= 130mW/7mr 2 . 



which gives a quadratic dependence with the field 
strength and is inversely proportional to the frequency. 
Note also that there is no dependence of the gap with 
the orientation of the k-vector. The polarization yields a 
maximum gap for the circularly polarized laser field while 
no net gap opens for the linearly polarized case. How- 
ever in Fig. [TJa noticeable corrections around the Dirac 
point are observed even for linear polarization. A deeper 
analysis shows that this is due to the field-induced lifting 
of the degeneracies between the states |k, ±1) ± around 
e — 0. This effect is of the same order as the correc- 
tion given by Eq. [3] and becomes stronger as the photon 
energy is reduced. 

Up to now, we have shown how the laser field modifies 
the electronic structure of 2d graphene that would be ob- 
served in any experiment carried out over a time much 
larger than the period T. A relevant question is if these 
effects would be observable in a transport experiment. To 
such end we compute the transport response at zero tem- 
perature using Floquet theory |26| applied to a 7r orbitals 
Hamiltonian. The electromagnetic field is accounted 
through the Pcicrls' substitution which introduces an ad- 
ditional phase in the hopping 7y connecting adjacent 

sites Ti and r,-: 7^ = 7oexp (i^ f* 3 A(t) • drj, where 

70 ~ 2.7 eV is the hopping amplitude at zero field and 
$0 is the quantum of magnetic flux. For computational 
convenience we use an armchair edge structure with the 
vector potential A(i) = A x cos fit x + A y sin fit y. 

Calculating the two-terminal dc component of the con- 
ductance as a function of Fermi energy is not, a priori, 
computationally easy. An efficient solution is achieved by 
decomposing the Hamiltonian into independent transver- 
sal channels as explained in . Although this decompo- 
sition is preserved by the interaction with the laser only 
for linear polarization (A y = 0), it is already enough to 
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give a flavor on the transport effects. Figure [3] gives the 
computed dc conductance 0, [26| for a stripe 1 //m wide 
and f /im long. One can appreciate in Fig. [3] that the 
strong depletion areas observed in the DOS of Fig. [1] 
are indeed mirrored in the dc component of the conduc- 
tance. Accordingly, the predicted dynamical gaps can be 
unveiled through transport measurements. 

In summary, we have shown that it is possible to use 
a laser field in the mid-infrared range to tune the elec- 
tronic structure of graphene and its electrical response. 



The modifications are predicted to arise both around the 
Dirac point and at ±M2/2. Moreover, since the results 
arc strongly dependent on the laser polarization, it may 
be used as a control parameter. We encourage experi- 
mentalists to pursue this exciting line of research. 
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